† Corresponding author. E-mail:
Project supported by the National Natural Science Foundation of China (Grant Nos. 11972177, 11972122, 11802103, 11772144, 11872030, and 11572034) and the Scientific Research Starting Foundation for Scholars with Doctoral Degree of Guangdong Medical University (Grant Nos. B2019042 and B2019021).
For conservative linear homogeneous nonholonomic systems, there exists a cotangent bundle with the symplectic structure dπμ ∧ dξμ, in which the motion equations of the system can be written into the form of the canonical equations by the set of quasi-coordinates πμ and quasi-momenta ξμ. The key to construct this cotangent bundle is to define a set of suitable quasi-coordinates πμ by a first-order linear mapping, so that the reduced configuration space of the system is a Riemann space with no torsion. The Hamilton–Jacobi method for linear homogeneous nonholonomic systems is studied as an application of the quasi-canonicalization. The Hamilton–Jacobi method can be applied not only to Chaplygin nonholonomic systems, but also to non-Chaplygin nonholonomic systems. Two examples are given to illustrate the effectiveness of the quasi-canonicalization and the Hamilton–Jacobi method.
The linear differential constrained system[1–3] is an important research object of analytical mechanics.[4–11] And linear differential constraints exist widely in many engineering problems, such as robots, bicycles, skateboards, electromechanical coupling systems, and so on. Many of these linear differential constraints are homogeneous.[12–24]
Linear differential constrained systems include holonomic systems and linear nonholonomic systems. When all differential constraints are integrable, the linear differential constrained system is called a holonomic system, and can be reduced into a geometric constrained system. When some differential constraints are not integrable, the linear differential constrained system is a nonholonomic system. In classical analytical mechanics, the correct motion equations of linear differential constrained systems can be written by the multiplier method. But because the motion equations with multipliers is a differential-algebraic equations, it is usually difficult to solve directly.
However the motion problem of conservative holonomic systems has been perfectly solved. In this process, the canonicalization of holonomic systems is a key. By introducing generalized coordinates qμ, the motion equations of conservative holonomic systems can be firstly reduced to the Lagrange equation without multiplier
![](cpb_29_6_064501/cpb_29_6_064501_ieqn1.gif)
From the geometric point of view, the canonicalization shows that a holonomic system has a natural symplectic structure.[27–34] If the symplectic form on the cotangent bundle of a holonomic system is recorded as Ω = dqμ ∧ dpμ, the Hamilton’s equations (
It is generally believed that nonholonomic systems have not the natural symplectic structures unless the Helmholtz condition is satisfied. For example, if the nonholonomic constraint of a dynamical system is
Because quasi-velocities is more consistent with geometric characteristics of the nonholonomic constraint than generalized velocities, applying nonholonomic frame, which consists of independent quasi-velocities
![](cpb_29_6_064501/cpb_29_6_064501_ieqn2.gif)
In this paper, by the first-order linear nonholonomic mapping theory,[40–43] we will give a method to realize the quasi-canonicalization of linear homogeneous nonholonomic systems. That is, by introducing a set of suitable quasi-coordinates πμ and quasi-momenta ξμ, the motion equations of conservative linear homogeneous nonholonomic systems can be written as
In this paper, we call the canonicalization of nonholonomic systems quasi-canonicalization. This is to distinguish it from the canonicalization of holonomic systems, and to emphasize that the canonicalization of nonholonomic systems depends on the introduction of suitable quasi-coordinates πμ.
In the third section of this paper, it can be seen that the application scope of the quasi-canonicalization method is linear homogeneous differential constrained systems. Therefore, the quasi-canonicalization is actually an extension of the canonicalization of holonomic system. That is, we extend the canonicalization from holonimic systems to linear homogeneous differential constrained systems which include holonimic systems and linear homogeneous nonholonomic system. It is easy to prove that for holonomic systems, quasi-coordinates πμ and the quasi-canonicalization can be reduced respectively to generalized coordinates qμ and the canonicalization mentioned in the third paragraph above.
Obviously, the quasi-canonicalization provides a new way to extend the Hamilton–Jacobi method and symplectic algorithms from holonomic systems to nonholonomic systems. Due to the symmetry of Chaplygin nonholonomic systems, the Hamilton–Jacobi method of Chaplygin nonholonomic systems has been proposed based on the symmetry reduction.[44–46] But the quasi-canonicalization method does not depend on the symmetry of nonholonomic systems. So the Hamilton–Jacobi method based on the quasi-canonicalization can be applied not only to Chaplygin systems, but also to non-Chaplygin systems. An example of studying non-Chaplygin systems by the Hamilton–Jacobi method is provided in the fifth section of this paper.
The structure of this paper is organized as follows: the first-order linear nonholonomic mapping theory is briefly introduced in the second section of this paper. The quasi-canonicalization method for linear homogeneous differential constrained systems is proposed in the third section. In this section, we firstly study the motion equations of linear homogeneous differential constrained systems on Riemann–Cartan configuration spaces. On this basis, we point out that there are nonholonomic mappings, by which we can write the motion equations of linear homogeneous differential constrained systems as Eqs. (
The Einstein’s summation convention is used throughout in this paper, such as i,j,k,l = 1,2,…,n; μ,ν,ρ,σ = 1,2,…,m; α,β = 1,2,…,n − m. The same kind of geometric quantities are represented by the same letter symbol, and the same kind of geometric quantities in different spaces are distinguished by different subscripts and superscripts. For example, gij and gμν represent respectively the metric of the space M and the space Π in this paper.
The idea of the first-order linear mapping method comes from the work of Kleinert et al.[47–51] In order to construct a Riemann–Cartan space with torsion, they proposed an embedding method by which a Riemann–Cartan space can be embeded into a known Euclidean space. Later, we extended the embedding method to the theory of first-order linear mappings.
Here we only briefly introduce the first-order linear mapping from flat spaces to Riemann–Cartan spaces. More about the first-order linear mapping theory can be seen in Refs. [52–55].
Take an n-dimensional flat space M with coordinates xi. The metrics gij and the connections
![](cpb_29_6_064501/cpb_29_6_064501_ieqn4.gif)
![](cpb_29_6_064501/cpb_29_6_064501_ieqn5.gif)
![](cpb_29_6_064501/cpb_29_6_064501_ieqn6.gif)
![](cpb_29_6_064501/cpb_29_6_064501_ieqn7.gif)
The space with a set of quasi-coordinates πμ is called the image space of the mapping (
It can be seen from Eq. (
When the torsion of Π is zero, Π will reduce to a Riemann space. In this case, Π is a submanifold of M.
It can be seen that Corollary
In this section, we study the motion equations of dynamical systems on Riemann–Cartan configuration spaces. On this basis, we can realize the quasi-canonicalization of linear homogeneous nonholonomic systems by the first-order linear nonholonomic mapping satisfying the condition (
First, let us briefly review some basic knowledge. When we study a dynamical system composed of N particles, the configuration space of the system is an (n = 3N)-dimensional flat space M with coordinates xi. the metrics of M are
If the system is subjected to n – m first-order linear homogeneous differential constraints
If an independent first-order linear differential equations are defined as
![](cpb_29_6_064501/cpb_29_6_064501_ieqn17.gif)
Next, we study the various motion equations of constrained systems in the reduced Riemann–Cartan configuration space Π. Especially, we will provide the quasi-canonical equations for conservative linear homogeneous nonholonomic systems in Proposition
The right side of Eq. (
According to Corollary
From the geometrical point of view, the reduced configuration space Q with generalized coordinates qμ always is a Riemann space. So there always is a natural cotangent bundle T*Q with canonical coordinates (qμ, pμ) for a holonomic system. This is the geometric basis of canonicalizations of holonomic systems. Obviously, the canonicalization process is independent of the selection of generalized coordinates qμ.
From Corollary
From the geometrical point of view, for a linear homogeneous nonholonomic system, the reduced configuration space Π is generally a Riemann–Cartan space with torsion. Because of the torsion, there is not a cotangent bundle on Π. Therefore, in order to achieve the quasi-canonicalization of the nonholonomic systems, a set of suitable quasi-coordinates must be defined, so that the reduced configuration space Π has no torsion. In this way, we can define the cotangent bundle T*Π with canonical coordinates (πμ,ξμ). This is the geometric basis of quasi-canonicalizations for linear homogeneous nonholonomic systems.
T*Π and T*Q are both submanifolds of the cotangent bundle T*M with canonical coordinates (xi,yi). The difference is that T*Q is embedded, while T*Π is only immersed. The symplectic form on T*Π (or T*Q) is the restriction of dxi ∧ dyi on T*Π (or T*Q). From
As an application of the above quasi-canonicalization, we study the Hamilton–Jacobi method for linear homogeneous nonholonomic systems.
Obviously, if the mapping (
For a linear homogeneous nonholonomic system, since the mapping (
In order to obtain the first integrals expressed by xi and
Of course, from the derivative of Eqs. (
![](cpb_29_6_064501/cpb_29_6_064501_ieqn23.gif)
We study the motion of the particle in the following two case:
The potential energy is U = x2. The potential energy is U = (1/3)(x2)3.
(I) When U = x2, the motion of the constrained particle can be solved by the traditional method. We first solve the problem by the multiplier method, then apply the quasi-canonicalization and the Hamilton–Jacobi method to solve it. It can be seen that the results of the two methods are equivalent.
The Lagrange equations with a multiplier λ of the particle are
Next, we apply the quasi-canonicalization method to solve this problem.
We introduce the following nonholonomic mapping:
The quasi-Hamilton’s equations of the system is
We calculate the derivative of Eqs. (
It can be seen that the solution (
(II) When
The motion equations of the particle with a multiplier is
Next, we try to solve this problem by the quasi-canonicalization and the Hamilton–Jacobi method. We also introduce the nonholonomic mapping (
The quasi-Hamilton’s equations is
In order to obtain the first integrals expressed by xi and
![](cpb_29_6_064501/cpb_29_6_064501_ieqn26.gif)
Equations (
The problem is completely solved by the Hamilton–Jacobi method.The correction of the solution (
Obviously, the constraint (
We introduce the following nonholonomic mapping:
The quasi-Hamilton’s equations is
In order to obtain the first integral expressed by xi and
![](cpb_29_6_064501/cpb_29_6_064501_ieqn29.gif)
The second equation of Eqs. (
The quasi-canonicalization of constrained systems depends on finding a set of suitable quasi-velocities. It is natural. The partial differential equations (
From the above two examples, we can see that
How to extend the quasi-canonicalization method in this paper to other differential constrained systems (such as linear nonhomogeneous differential constrained systems) is also worthy of further study. For linear nonhomogeneous differential constrained systems, the difficulty of the generalization lies in that the image space Π of the first-order linear mapping implying constraints has a time dimension. When the time dimension and the space dimension cannot be separated, the constraint force cannot be geometrized.[43]
We extend the canonicalization and the Hamilton–Jacobi method from holonimic systems to linear homogeneous nonholonomic systems in this paper. It means that the motion equations of conservative linear homogeneous nonholonomic systems can be written into the canonical form with quasi-coordinates πμ and quasi-momenta ξμ. The canonicalization realized by quasi-coordinates πμ is called the quasi-canonicalization in this paper. In order to achieve the quasi-canonicalization, a set of suitable quasi-coordinates πμ need to be defined by a first-order linear mapping, which implies the constraints, so that the reduced configuration space Π of linear homogeneous nonholonomic system is a Riemann space with quasi-coordinates πμ. Since there is naturally a cotangent bundle T*Π with the symplectic form d πμ ∧ dξμ on the Riemann space Π, the motion equations can be written as the canonical form by the set of suitable quasi-coordinates πμ and quasi-momenta ξμ. And the Hamilton–Jacobi method for linear homogeneous nonholonomic systems is a natural result of the quasi-canonicalization. Because the Hamilton–Jacobi method based on the quasi-canonicalization does not depend on the symmetry of linear homogeneous nonholonomic systems, it is suitable for non-Chaplygin nonholonomic systems. How to extend the quasi-canonicalization and the Hamilton–Jacobi method to more general nonholonomic systems is worthy of further study.
[1] | |
[2] | |
[3] | |
[4] | |
[5] | |
[6] | |
[7] | |
[8] | |
[9] | |
[10] | |
[11] | |
[12] | |
[13] | |
[14] | |
[15] | |
[16] | |
[17] | |
[18] | |
[19] | |
[20] | |
[21] | |
[22] | |
[23] | |
[24] | |
[25] | |
[26] | |
[27] | |
[28] | |
[29] | |
[30] | |
[31] | |
[32] | |
[33] | |
[34] | |
[35] | |
[36] | |
[37] | |
[38] | |
[39] | |
[40] | |
[41] | |
[42] | |
[43] | |
[44] | |
[45] | |
[46] | |
[47] | |
[48] | |
[49] | |
[50] | |
[51] | |
[52] | |
[53] | |
[54] | |
[55] |